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Abstract 



We asymptotically analyse the volume-random variables of general, symmetric 

and cyclically symmetric plane partitions fitting inside a box. We consider the re- 

CN spective symmetry class equipped with the uniform distribution. We also prove area 

limit laws for two ensembles of Ferrers diagrams. Most of the limit laws are Gaussian. 



^ 1 Introduction 

^ A plane partition fitting inside an (r, s,t)-box is an r x s-array 11 of non-negative integers 

O Pi,j ^ t with weakly decreasing rows and columns. It can be visualised as a pile of unit 

j> cubes in the box B{r, s, t) := [0, r] x [0, s] x [0, t] "flushed into the corner", see flgure[l]and 

k>J [2] ■ The flgure also illustrates the connection to tilings of a hexagon of side lengths r, s and 

^ t by lozenges. Yet another interpretation is viewing such a pile of cubes as an order ideal 

in the product of three flnite chains (total orders) with the respective lengths r, s and t. 

In the following we mean by "plane partition" one which flts inside an (r, s, t)-box. The 

volume of a plane partition is the sum of its parts or the number of unit cubes in the pile 

or the cardinality of the order ideal, respectively. If one side length, say t of the bounding 

box is one, such a pile can be regarded as Ferrers diagram fltting inside an r x s-rectangle 

of area equal to the volume. Let VV{r,s,t) denote the set of all plane partitions fltting 

inside B{r, s,t). A plane partition 11 G VV{r, r, t) is called symmetric if the corresponding 

pile of cubes is symmetric about x = y. To put it differently the cube {i,j, k) belongs to 

the pile, if and only if {j,i,k) does. The respective tiling is symmetric w.r.t. a vertical 

axis. A plane partition 11 G VV{r, r, r) is called cyclically symmetric if the corresponding 

tiling is invariant under a rotation about 27r/3, i.e. the cube {i,j, k) belongs to the pile if 

and only if (fc, i,j) and (j, k, i) do. Call these subsets SW{r, t) and CSW{r) respectively. 
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Figure 1: A plane partition in (4, 5,4)-box 



Several authors discussed properties of randomly chosen plane partitions, for example in 
[9] diagonal sums are considered. In [?] an "Arctic Circle Theorem" for tilings of a large 
hexagon is proved, which states that a typical tiling looks periodic near the corners of the 
hexagon and unordered inside the inscribed ellipse. 

It is natural to ask also for volume limit laws as the generating functions counting plane 
partitions by volume are available for the above classes. For the following formulae we 
refer to [2]. The volume generating function of VV{r,s,t) is 
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Symmetric plane partitions in B{r,r,t) have the volume generating function 
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Here the first product corresponds to singleton orbits and the second product corresponds 
to doubleton orbits of B{r, r, t) under switching the first two coordinates. For the cyclically 
symmetric plane partitions we have the volume generating function 
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where the first product runs over C3 orbits of B{r,r,r) with all coordinates distinct, the 
second runs over orbits with exactly two coordinates equal and the third corresponds to 



singleton orbits. The above generating functions are polynomials in q whose degree is the 
volume of the bounding box. 

We equip VV{r,s,t), SVV{r,t) and CSWij-) with the uniform distribution ("uniform 
fixed bounding box ensemble"). By X = Xj-st (resp. = X^, Xr) denote the corresponding 
volume random variables. In the subsequent discussion we will drop subscripts as the side 
lengths are always r, s and t (resp. r and t, r). The probability generating function of this 
random variable is P{q) '■= G{q)/G{l) where G{q) is the respective generating function. 
In the VV{r, s, t) case this reads with the notion Uijk := i + j + k — 1 

In the following g{x)^^^ denotes the iV*'* derivative of g w.r.t. x. Moments of X of arbitrary 
order A^ exist and can be computed via 



E(X^) = (-/)^(P(e^^))^^^ , (1.5) 

x=0 

where P(e^^) (/^ = —1) is the characteristic function of X. 

2 Mean, variance and concentration properties 

Denote by /i, /i^pp and ficspp the mean value and by a, a^pp and a^spp the standard deviation 
of the volume variables on VV{r, s,t), SVV{r,t) and CSVV{r), respectively. 

Lemma 2.1. The volume distributions are symmetric about half the volume of the bounding 
box, and hence the expected volumes are ii = rst/2, fj,spp = rH/2 and ficspp = t'^/2. 

Proof. Let v be the volume of the bounding box B. We show that X and v — X are equally 
distributed. Let a plane partition of volume k be represented by a pile of k green cubes 
inside B. Now fill up the empty space in B with red cubes. The red cubes represent a plane 
partition of volume v — k. This construction describes a bijection between plane partitions 
of volume k and v ~ k which respects symmetry and cyclic symmetry. Now the lemma 
follows as E(X) =¥.{v-X). U 

In order to prove the Gaussian limits we consider the characteristic functions P(e^^) 
of the random variables X. More precisely, since sums are easier to handle than products, 
we will study the logarithms of the characteristic functions. The following lemma enables 
us to compute explicit formulas for the variances of the random variables and to estimate 
the Taylor coefficients of the functions in question. 

Lemma 2.2. For positive real numbers a, c with a > c>l, we have the expansion 

log f ^"7/^^'rf? l = E H,c{a)x-. (2.1) 

V a(l-e("-^)^) / ^^ '^ ' ^ ' 



Here Hn,c is a polynomial of degree N ~ 1 in a. In particular we have Hi^c{oi) = \ and 

2 

H2,c{o) = —^a + 1^. Furthermore there is a positive constant D, such that the inequality 

\HnM)\ < D ■ a""-' ■ {2cf (2.2) 

holds for all N eN. 
Proof. Define the function g{t) by 



g{t) = log 
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(2.3) 



Tlien tlie llis of (2.1) is easily seen to be equal to g (ax) — g {{a — c)x) . We have the 
following series expansion for g(t) : 



9{t) = log 
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Observe that the singularities of g of smallest modulus are ±27r/, so the numbers \bi\i\ 
decay like (27r)~^ and hence are bounded by some constant D. The iV*'* coefficient in the 
Taylor expansion of g (ax) — g {{a — c)x) about x = is in fact a polynomial of degree 
N — 1, namely 
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As a > c > 1 it can be estimated as follows: 
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This finishes the proof of Lemma 2.2 
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Now we can easily compute the variances. The formula for cr^ already appeared in [9]. 
Lemma 2.3. We have 

1. (7^ = j^(r^st + rsH + rst^) = ^rst{r + s + t), 



2. al^^ = |tr3 + It'^r^ - ^t^r + ^tr^ - \tr and 

^ ■ ^ cspp 4 ' 2 



Proof. Recall that the variance of a random variable Y can be obtained as the Y{Y) 
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— log (P(e'^^)) , where P{q) is the probability generating function of Y. In the VV{r, s, t) 

x=0 

case we apply this to P{q) as in (1.4) and obtain 

-iog(p(e-))"| =-EEE'-^*"'" 

i=\ j=i k=i 
where aijk = i+j + k — 1. According to Lemma 2.2, each summand on the rhs is equal to 

2\H2,iiaijk) ■ 
and a straightforward calculation yields 

"• " ^ /I 1 \ 1 

J2J2J2[q^'+^ + ^-'^^-Y2) =^{r'st + rsH + rse). (2.9) 
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The other variances are calculated in an analogous way with suitable choices of a and 

c. n 

Now we can investigate concentration properties of the families of volume random 
variables, i.e. we study the quotients of standard deviation and mean when the box gets 



large. According to Lemma 2A_ and Lemma 2^ these quotients tend to zero if at least two 
side lengths of the bounding box tend to infinity. The latter is in particular satisfied when 
r — > oo in the SVV{r, t) and CSVV{r) case. In the general case of VV{r, s, t) we have 
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If say r and s are unbounded, the right hand side of (2.10) tends to zero for r, s — ^ oo and 



the family is concentrated to the mean, i.e. for every S > we have 
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On the other hand, if two coordinates are fixed, say r, s, then the quotient is bounded away 
from zero by 5^ for t — > 00. These families are not concentrated. Similar considerations 
hold in the case of SVV{r, t) when r is fixed and t — i> 00. In the next section we investigate 
the limit laws in these two cases. 



3 Limit laws 

We first give a result for the concentrated families. 



Proposition 3.1. If at least two ofr,s,t tend to infinity, the family 

X- fi' 



Y :-- 
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of normalised volume random variables on VV{r,s,t) converges in distribution to a stan- 
dard normal distributed random variable. The same statement holds for the volume vari- 
ables on SVV{r,t) if at least r — > oo and on CSW{r) z/r — > oo. 

Proof. For VV{r,s,t). Let r, s — > cxd and t vary arbitrarily. Tlie cliaracteristic function of 
Y is (p{x) := e~^^^/'^P (e^^f'^^ , with P as in eq. (1.4). We prove that log (0(a;)) — > — ^ for 
X G M if r, s -^ oo. Then by Levy's continuity theorem |1] the assertion follows. 



By Lemma 2.1 Y and —Y have the same distribution. So the characteristic function is real 



valued and the coefficients for odd A^ vanish. By Lemma 2^ we have the Taylor expansion 
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(3.1) 



According to the estimate (2.2) we can bound the modulus of the 2A^*'^ summand, N > 2 
by 
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Plugging the explicit expression for o"^ of Lemma 2.3 into the rhs of (3.2) we obtain the 
estimate for the 2A^*'' summand, A^ > 2, in the expansion (3.1) 



j=l j=l fc=l 



xl 

a 



2N 



<D|x|2^-48 



iN 



1 




1 




1 


rg 


+ 


rt 


+ 


St 



7V-1 



(3.3) 



Summing in (3.3) over N >2 then yields the estimate 
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(3.4) 



The rhs tends to zero for any fixed real x if r and s tend to infinity. This proves the 
assertion for the VV{r, s, t) case. The other cases are shown with similar estimates. D 

Now we consider the non-concentrated case. If r = s = 1 we have a single column 
of unit cubes of height t. The volume clearly is uniformly distributed. So, if rs > 1 the 
volume random variable can be viewed as sum of dependent random variables with values 
in {0, . . . ,t}. Now the easiest guess is that for large t the dependence vanishes and the 
volume of a single column is uniformly distributed. This guess is the right one as the 
following proposition shows. 



Proposition 3.2. If r and s are fixed and t tends to infinity, the family 

Xrst 

~T 

of rescaled random variables converges in distribution to the {rs)-fold convolution of the 
uniform distribution on [0,1]. In the SVV{r,t) case the so rescaled sequence converges 
in distribution to the convolution of r factors of the uniform distribution on [0, 1] and 
r{r — l)/2 factors of the uniform distribution in [0, 2]. 

Proof. For VV{r,s,t). We show that the Fourier transform of Zt converges pointwise to 
{ ^ — j , which is the Fourier transform of the rs-fold convolution of the uniform distri- 



bution on [0, 1]. The Fourier transform of Zt is P{e^^^*), with P as in (1.4). We expand 
the product running from 1 to t in P{e^^^^). All but the last term in the numerator and 
the first term in the denominator cancel out: 

The single factors are easily seen to converge to ^ j~^ . 

The SVV{r,t) case is worked out analogously. D 

4 Ferrers Diagrams 

A Ferrers diagram is a convex lattice polygon which contains both upper corners and the 
lower left corner of its smallest bounding rectangle. When t = 1 a boxed plane partition 



can be viewed as a Ferrers diagram fitting inside an r x s rectangle. Then formula (1.1) 
reduces to the g-binomial coefficient ['^^**] . The class of Ferrers diagrams with h rows and 
w columns has the area generating polynomial 

Q'^'-'n-i'], (4.1) 

since at least the h + w — 1 unit squares which constitute the top row and the leftmost 



column are contained in such a polygon. So, by Proposition 3J^ a Gaussian area limit law 
arises in this uniform fixed height and width ensemble, when both height and width tend 
to infinity. 

In statistical physics [7] as well as in combinatorics pp, one is also interested in area limit 
laws in the so-called uniform fixed-perimeter ensembles, where all Ferrers diagrams of a 
fixed half-perimeter are considered equally likely. For fixed m > the area generating 



function of Ferrers diagrams of half-perimeter m + 2 is obtained by summing formula (4.1 ) 
over all pairs [h, w) G N^ with h + w = m + 2. This can be written as 

m r- 



h=0 



Observe that the index of summation in (4.2) can be interpreted as the height minus one 
and that the total number of such Ferrers Diagrams is 2™. So the following probability 
generating function (PGF) 
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describes the ensemble of Ferrers diagrams of half-perimeter m + 2 counted by height and 
area. The additional height parameter allows us to use the results for height- and width- 
ensembles dicussed above. More precisely, we condition the area variable in the perimeter- 
ensemble on the height variable. In what follows we will prove that the joint distribution of 
the (properly rescaled) height and area variables in this ensemble converge in distribution 
to the two-dimensional standard normal distribution. The following proposition from [0] 
(see also [8J) is taylor made for this situation. 

Proposition 4.1. Let Xm, Ym be real valued random variables and let Ym be supported in a 
lattice Lm '■= {am + kSm\k G Z} , where 6m > and am G K, i-e. f(Ym G Lm) = 1- Suppose 
Yjn satisfies a local limit law fi with a density g{y) w.r.t. the Lebesgue measure on M. (this 
implies 6m — > 0), i.e. for all y eM. and every sequence (y^) with y^ G L^ and ym —^ y we 
have F{Ym = ym)/6m -^ 9{y)- Suppose further that for fi-almost all y E M. the conditional 
distributions f{Xm G -{Ym = ym) converge weakly to a measure i^{-,y) (ym ^ y, Vm & Lm). 
Then the joint distribution of {Xm, Ym) converges weakly to the measure z/ defined by 
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for all Borel sets A, 5 C M. 

For q= 1 formula (4.3) is simply 2~'"m(1 + m)™, the PGF of the height random variable 



Hm- It satisfies the assumptions of |5i Theorem IX. 14] and hence a Gaussian local limit 
law arises. The mean and standard deviation of the height are easily computed to be 
asymptotically equal to m/2 and y/m/2, respectively. So let Ym = ^{Hm — m/2)/y/rn. 
The random variable Ym is supported in %j ^fm. Let y G M and ym -^ V-, where each 
ym £ 1>l ^fm. The area random variable Am, conditioned on the event \Ym = ym} , has 
PGF 
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An application of Lemma 2.1 and Lemma 2.3 with r = m/2+ym\/rn/2., s = m — r and t = 1 



shows that the mean and the variance of the conditioned area variables are asymptotically 
equal to m^/8 and m^/48, respectively. So according to Proposition 3.1, the rescaled area 
variable Xm '■= (Am — m'^/8)/^/m^^j48 conditioned to {Ym = Vm} converges weakly to a 



Gaussian random variable with mean zero and variance one. Now Proposition 4.1 yields 



Proposition 4.2. Denote Am and Hm the random variables of area and height of a Ferrers 



diagram in the uniform fixed perimeter ensemble described by (4.3). Let 

- m^/8 ^^ Hm - m/2 



X„ 
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Then, as m ^ oo, {Xm, Ym) converges weakly to the two-dimensional standard normal- 
distribution, n 

The last result can also be obtained by analysing a g-difference equation satisfied by 
the generating function of Ferrers diagrams counted by height, area and half-perimeter. 
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